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1. Introduction 

It has been customary to accept as the basis for a seasonal 
adjustment procedure a model for the seasonal component 
of a time series which has the form: 

n 

s t =Y,^\ (i.i) 

t=i 

where 

St = seasonal component at time t 
a.j = complex-valued coefficients 
ojj = seasonal frequencies such that u>j = jto i 
to i = l cycle per year. 

n = a constant, such that there are 2 n equally spaced 
observations per year.* 

This is the usual finite Fourier series decomposition of an 
arbitrary shaped periodic function of period 

UJi 

The model for the observed time series is then given 
as: 

X t = S t + Z u (1.2) 

where 

X t = observed time series 

Z t = a periodic stochastic time series. 

It frequently appears to be appropriate to consider the 
logarithms of the original variables when postulating an 
additive model of this form. However, we will not discuss 
this issue here. 

The estimation problem for this process is that of es¬ 
timating S t . from a realization of X t . The estimation of 
St. permits the computation of Z t , which is generally sup¬ 
posed to reveal certain information which is not evident 
in X t . For the case where Z t is a wide sense stationary 
process and the a/s of equation (1.1) are constants in¬ 
dependent of time, the estimation problem seems to be 
completely solved (Hannan (1963), Jorgenson (1964)). 

However, statisticians (Wald (1936), Shiskin (I960)) 
have long recognized that observed “seasonal” patterns in 
economic time series do not conform to the fixed, strictly 


periodic, form given by equation (1.1). This observation 
has led to attempts to modify this basic model. These 
attempts have followed two principal lines. Both lines in¬ 
volve making the aj's of equation (1.1) functions which 
change with time. 

The first line of development has resulted in the pos¬ 
tulate that each ay is a function of time only. Estima¬ 
tion requirements obviously demand that each ay be a 
slowly varying function time when compared to cos(u>i). 
The estimation procedures for this model simply involve 
estimating the aj's as constants over short successive time 
intervals and then constructing the time series from 
these estimates. The pioneering practical work of Shiskin 
(1960) and his co-workers is well known in this area as is 
the analytic work of Hannan (1964). 

An interesting point to observe about this method is 
that it nearly preserves the basic linearity of the process 
which is implied in equation (1.1). In the case where the 
a/s are allowed to vary slowly with time, but indepen¬ 
dently of each other, we are assuming that the seasonal 
is still generated by a linear, but mildly non-stationary 
process. 

The other line of development has involved making 
the aj's both functions of time and of each other. Only 
the average values of the a/s over the realization are es¬ 
timated using equation (1.1). In the procedure as orig¬ 
inally developed by A. Wald (1938) the average values 
were constants, however (Godfrey and Karreman (1964)) 
the average values may also be estimated as independent 
functions of time in order to deal with certain kinds of 
non-stationarity. After the average values of the aj's have 
been estimated it is assumed that over relatively short 
time-intervals the aj's vary proportionately. Thus the ra¬ 
tios of the aj's remain nearly constant, but the aj's vary 
to take account of the changing amplitude of the seasonal 
variation. 

This second line of development takes account of a 
certain kind of non-linearity by removing the assumption 
of independence when estimating the variation of the aj's. 
However, the method still starts with the basic linear mo¬ 
del in order to estimate the mean values of the aj's, and 
then it imposes a rather special character on the variation 
of the aj's. 

The seasonal adjustment method to be proposed here 


* 2 n is used because the highest frequency which can be 
discriminated given an equispaced record at time points 
t = 1,2,..., m is lo„ = 7r. 
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is developed on the basis of a rather more general model 
concerning the process by which the seasonal component 
is generated. The model will be shown to include the two 
models indicated above as special cases. 


2. The Model of Seasonal Variation 

The basic assumption of the procedure to be described 
below is that the seasonal variation is produced by a few 
unobservable factors. These factors may be thought of as 
time series which determine the time paths of the time- 
dependent a/s of ecfuation (1.1). Thus the estimation 
problem is that of estimating these unobservable factors. 
The estimation of the factors may follow a time series gen¬ 
eralization of the classical principal factor method of factor 
analysis. The number of (orthogonal) factors required will 
depend on the explanatory power of each additional fac¬ 
tor. Clearly, the power of the method will depend on how 
many factors are required to explain the seasonal varia¬ 
tion. If only one factor is required, then the method will 
form a a simple generalization of Wald’s method (Wald 
(1936)). If n factors are required, then the method will be 
very similar to Hannan’s method (Hannan (1963)) as the 
a/s will be independently estimated. 

We may now proceed to the formal description of the 
model and the technique of its estimation. We assume the 
usual additive model for the observed time series 

X t = S t + Z t (2.1) 

where 

X t = observed time series 
St. = seasonal variation 

Z t . = time series with continuous spectrum whose den¬ 
sity is smooth in the neighborhood of season frequen¬ 
cies. 

Thinking in terms of the frequency decomposition of 
St, we nay represent S t by: 

n 

s t = Y, x ^ eiUit ( 2 - 2 ) 

i—1 

where the complex-valued time series x^t each have spec¬ 
tral power only in the neighborhood of zero frequency. We 
now postulate that the x itt are not independent series, but 
rather they are generated by a smaller set of p (p < n) se¬ 
ries. Therefore, we have: 

p m 

A k,i{j)yk,t—j ( 2 - 3 ) 

j =0 

where 


Dk,t-j = form, for l = l,...,pa set of time series, 
m = maximum number of lags used, 

Ak,i(j) = matrix of coefficients. 

Each series x^t is generated by a linear combination 
of past values of some or all of the y k ,t- 

Instead of estimating the Xij as if they were indepen¬ 
dent, it would appear reasonable to attempt to estimate 
the smaller set of variables y k ,t- We may therefore write 
equation (2.3) as 

p m 

Xi = ^2 ^2 A kA z 3 y k (2.4) 

k= 1 j— 0 

using the z-transform operator defined by: 

OO 

x = 

t=0 

We then represent the quantities A kt iZ 3 by B k ,i and 
solve the equation (2.4) for the y k 's. Thus we have 

p 

Vk = £ C k ,iXi (2.5) 

2—1 

where 

Ck,i = B ki 

We may then consider the estimation of y k by first com¬ 
puting y k , t by: 

p m 

Vk,t — ^ ^ ) Ck : i{j)Xi : t.—j T Ciy(/l) (2.6) 

i— 1 j= 0 

where: £i,t(/c) = independent random numbers. 

We then estimate by 

Vi,t = D{ili,t) (2.7) 

where: D{ ) is a low pass filter. 

While equation (2.3) states that the n series of sea¬ 
sonal coefficients x- kt are determined by a set of p time 
series yt^t, we will want to estimate the y^t from equations 
(2.6 and 2.7). We will then use an equation of the form of 
equation (2.3) to get estimates of the x ijt . This equation 
may be stated as: 

p m 

Ak,i(j)yk,t-j ( 2 . 8 ) 

k—lj—0 

We will now describe the estimation procedure in 
somewhat greater detail. 
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3. Estimation of the Model 

The estimation procedure to be presented here is one which 
seams natural, but it is only a first attempt at estimation 
and by no means represents a complete solution of the 
estimation problem. Further work, both theoretical and 
practical, is required in this general area. 

The first step in the estimation procedure is to com¬ 
pute the time series x by complex demodulation. The 
Xi t t are thus given by: 

Xi, t = L(x t e lu>it ) (3.1) 

where 

L( ) - denotes a low-pass filter with bandwidth of 

about wi. 

The Xi t t form estimates of the complex-valued Fourier co¬ 
efficients at frequency w*. In effect, each term in the series 
Xj i t is an estimate over the time interval t — m\ to t — m 2 
where m = mi +m. 2 , and m is the number of terms in the 
low-pass filter. 

These x lt may be considered as linear estimates of 
the seasonal variation. If we set m equal to the number 
of observations (and mi = m 2 = m/2) we then have a 
single (complex-valued) number for each series Xi t t- These 
are clearly estimates of the Fourier coefficients in the sta¬ 
tionary linear model. Letting m be less than n gives us 
estimates of time varying Fourier coefficients which are 
similar to what the Census Method [Shiskin (1960) and 
see also the relevant analysis of Rosenblatt (1963)] and 
the BLS method produce. As m is decreased, for a given 
n, the variability of x-^t increases. Depending on one’s as¬ 
sumptions about the variability of the amplitudes of the 
seasonal there will be some m which minimizes the ex¬ 
pected error of the estimate over the realization of n ob¬ 
servations. 

In forming estimates of time varying coefficients it has 
been general practice to make mi = m 2 - This implies an 
indifference toward the direction of time. For the purpose 
of forecasting future values of Z t it would seem more natu¬ 
ral to set m 2 < 0. In this case we are dealing with a linear 
prediction problem and the “best” operator may be the 
Wiener optimum predicting filter. 

However, for our present purposes the appropriate fil¬ 
ter is simply one which removes the frequency terms above 
cos(wif) which will appear in the Xj, it e lUit series before fil¬ 
tering. For this purpose a low order (say 8 or 10 term) 
weighted moving average should be adequate. 

Given these complex-valued series we now wish to 
estimate the transformation which will give the p series 
which explain most of the variation in the n series. A nat¬ 
ural way to do this is to estimate the principal factors at 


a set of modulating frequencies. We may then combine 
these factors in order to form the series iji,t as given in 
equation (2.6). 

The modulating frequencies which should be chosen 
for this estimation are naturally A = 0, n/kn, 2n/kn ,..., 
7T /n where k may equal, say, 2 or 3. 

The number of modulating frequencies which we can 
hope to estimate is severely limited by the amount of data 
available. The frequency n/kn has a period of 2k years. 
Thus, for k = 3 we are attempting to estimate a mod¬ 
ulating frequency with a period of 6 years. In order to 
estimate the principal factors we band-pass filter the se¬ 
ries Xi t t around A j to obtain a series Xi t t(Xj). We then 
compute the covariance matrix: 



where: 

1 N 

x i,q{ X j) = ^^2 X iA X j) X *q,t( X ^ ( 3 ' 2 ) 

v t -1 

The eigenvectors associated with the p largest eigenval¬ 
ues of R(Xj) define the linear combination of the p series 
yij which best explain the variation of the Xi,t at the fre¬ 
quency A j. We therefore have a p x n matrix, F(Xj), at 
each frequency A j. 

From these matrices we can compute the coefficients 
C k ,iU) °f equation (2.6). These coefficients will then give 
the linear combination of the x' i t s which explain most of 
the variance of the factors 

We can now derive estimates of the variation in the 
seasonal factors, iji.t, on the basis of assumptions concern¬ 
ing the model which generates these factors. The most 
natural assumption is that the are generated by a low 
order autoregressive model. Thus, we estimate the sea¬ 
sonal factors by: 

1 

Vi,t = £ K s )Vi,t-s (3.3). 

S=1 

The value of l and the values of the autoregressive coef¬ 
ficients b(s) should be determined after inspection of the 
yi t t and, probably, computation of their spectra. The “op¬ 
timum” values of the b(s) depend on the exact assump¬ 
tion concerning the y,^ and the amount of data available. 
However, in practice a standard set of b(s)'s will almost 
certainly give nearly as good results. 

We can now compute the according to: 

p m 

k—1 j =0 


(3.4) 
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where the Ak,i(j) is derived from the Ck,i{j ). 

The estimated seasonal series is then given by rernod- 
ulating and summing the x.^t : 

n 

= (3.5) 

i—1 

Finally, the seasonally adjusted series is given by: 

Z t = X t - St. (3.6) 


4. Relation of the Model to Linear processes 
and to Wald’s Model 

In order to explore the nature of this model somewhat fur¬ 
ther it is interesting to consider the expected result if the 
procedure is applied when in fact the seasonal was gener¬ 
ated by a linear process. It will also be interesting to see 
that special assumptions concerning the form of the matrix 
Ak : i(j ) or equivalently F(Xj) lead to a procedure which is 
similar to the Wald method of seasonal adjustment. 

4.1 Linear Seasonal Variation 

We may consider the effect of estimating by this proce¬ 
dure if the series S t is composed of n independent series 
with power at frequencies u>i for i = 1,..., n such that the 
Fourier coefficients vary independently with time. In this 
case it is to be expected that n factors would be required 
to explain the variation of the n series Xj,t and that these 
factors could simply be the series x^t- Thus, we are led, in 
the linear case, to estimation of the time varying Fourier 
coefficients by n independent autoregressive schemes. 

4.2 Relation to Wald’s Method 

We may also consider the case which, in a sense, is at the 
other extreme from the linear case. This is the case in 
which the Xj,t are given by: 

Xi tt = a,iX i, t i = 2,... ,n (4-2.1) 

This restriction on the x, h t implies that they are perfectly 
correlated and vary proportionally without any lag. Wald’s 
assumptions concerning the “pattern” of seasonal varia¬ 
tion imply this assumption. (Wald (1936) pg. 78; and 
Godfrey and Karreman (1964) pg. 15.) If equation (4.2.1) 
holds then we would expect that we could explain the vari¬ 
ation in the series by a single factor y for which the 


matrix Aky(j) would be given by: 

A k,i(j) = [l, a 2,---,a n ] j = 0, k = 1, ( 42 2) 

= o, j r 0 

If we were to estimate y by a 12 month moving 
average we would be exactly paralleling Wald’s method. 


5. Conclusion 

The method proposed here has been motivated by the ob¬ 
served behavior of a considerable number of economic time 
series. The most prominent characteristic of many eco¬ 
nomic time series, as originally observed by Wald (1936), 
is that the seasonal pattern seems to be stable but am¬ 
plitude modulated. This results in the amplitude of the 
fundamental and all the the harmonics varying together. 
Quite naturally the method presented here is similar to 
Wald’s method and, in fact, it includes his method as a 
special case. However, as the method is somewhat more 
general than Wald’s method it will provide reasonable es¬ 
timates in situations such that the modulation of the pat¬ 
tern is not present, but the amplitude of each harmonic 
varies independently. 

This method is far from completely developed or de¬ 
scribed in this paper. The next natural step is to program 
the method and study its behavior along the lines of the 
analysis in Godfrey and Karreman (1964). 
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